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MOMENTS OF INCLINATION ERROR DISTRIBUTION 
COMPUTER PROGRAM 


By T. R. Myler 
Vought Corporation 


SUMMARY 

This report descrihes a FORTRAN coded computer program which calculates 
orbital inclination error statistics using a closed-form solution. This 
solution uses a data base of trajectory errors from actual flights to predict 
the orbital inclination error statistics. The methods used to generate the 
error statistics are of general interest since they have other applications. 
Included in this report are program theory, user instructions, output 
definitions, subroutine descriptions and detailed FORTRAN coding information. 
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1.0 INTRODUCTION 


Since 1969f flight experience has been used as the basis for predicting 
Scout orbital accuracy. The data base of flight experience consists of errors 
in the trajectory parameters at orbit insertion as observed on Scout flights. 
This data base is used in a Monte Carlo analysis to calculate error statistics 
in the orbital parameters. The need arose to generate orbital inclination 
error statistics using a closed-form solution which included the option of 
considering non-zero mean values and correlation between the variables in the 
data base. To meet this need, the classical approach of expanding the orbital 
inclination relationship by a multivariable Taylor series was chosen. The 
mechanization of this process resulted in computer program Moments of 
Inclination Error Distribution (acronym MIED) and is described herein. 
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2.0 DEFINITIONS 


2.1 Notation 


Symbols used in this report are listed below with their definition 
and units. 

English Alphabet 


Ce 

east component of position, ft 

On 

north component of position, ft 


crossrange, n.mi. 

E(x) 

expected value of x 

1 

orbit inclination, deg. 

n 

number of samples 

R 

range, n.mi. 

^e 

earth radius, ft 

Greek Alphabet 

a 

longitude, deg 

A 

deviation from nominal 

f 

azimuth, deg. 

X 

latitude, deg. 

At 

statistical moment 

a 

standard deviation 


Superscripts 

(•) 

Other 

d X 
8 y 


adjusted value 


partial derivative of x to y 


- 3 - 



2 . 2 Flowchart Conventions 


Flowchart conventions used in this report are as follows 


Process 



Input/Output 



Subroutine 


Decision 




3.0 PROGRAM DESCRIPTION 


This section describes the utilization of the Scout data base, program 
theory, user instructions and output definitions. 


3.1 General 


The purpose of computer program MIED is to calculate orbital 
inclination error statistics from a data base of errors resulting from actual 
flights. The Scout flight history data base consists of orbit insertion 
errors in the trajectory parameters - altitude, velocity, flight path angle, 
flight azimuth, latitude and longitude. This data base is shown in Table 3.1* 

The Scout data base includes flights with launch azimuths ranging from 
easterly to southerly to slightly west of south. Latitude errors, which 
produce inclination errors, result from crossrange errors on easterly flights 
and range errors on southerly flights. Therefore, latitude errors of the 
flight history data base are not a consistent sample set of errors for a given 
flight azimuth at insertion. Since it is necessary to have a consistent set 
of sample latitude errors, the data base latitude errors are adjusted for the 
flight azimuth of interest. Knowing the nominal flight azimuth, latitude and 
their errors on each flight of the data base, range and crossrange errors are 
calculated. Since range and crossrange errors are independent of the flight 
azimuth, these errors can be converted to latitude errors for the flight 
azimuth of interest. This technique is detailed below: 



AR = 
ACR 


r^AX 

Fg Aa cos X 
Cjyj cos f + sm f 
= sin f — Cg cos 


f 


The above relationships are evaluated for each flight sample of the data base, 
resulting in range and crossrange errors for each flight. Note that latitude 
and longitude errors are needed to calculate range and crossrange errors. 

Since the desired end result is inclination error statistics for a specific 
flight azimuth and latitude at insertion, the above range and crossrange 
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errors are converted to a latitude error applicable to this specific flight 
azimuth, as shown below. 


C'|vj = ARcosf'+ ACR sin f' 

AX' = C 

where the primed values pertain to the conditions at which the 
inclination error statistics are desired. 


The above process yields a sample set of adjusted latitude errors and are used 
in conjunction with the azimuth errors of each flight to generate latitude and 
azimuth error statistics as described in Section 3-2. 


3-2 Program Theory 


The classical approach to generating statistical moments of a 
multivariable system is to use a Taylor series expansion. The derivation of 
the resulting relationships for a system of correlated variables is presented 
in Appendix 7B of Reference (l). Applying the orbital inclination 
relationship to the equations of Appendix 7B yields the algorithm used in 
computer program MIED. 

Orbit inclination has the following relationship to trajectory 
parameters : 


cos I = cos X sin f 


Thus, latitude and azimuth are the system variables. This equation is applied 
at orbit insertion to obtain the final orbit inclination. Adapting the above 
relationship to the equations given in Appendix 7B for two variables yields 
the following: 


Expected value = E(i) = cos jcos [ E(X)] sm [E(f)]| 

02, o 92, 


+ 'h 


a\2 

32 | 

axar 


+ 'h- 


|[X- E (X)l [r- E (f)]} 
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2 


2 

Variance - 


Third moment - 


Fourth moment 




9i 3i 


2 

3X 

ar 

3i 

a2, 

1a 

ax2 

3i 

a2, 

r* 

IT 

9f2 

32 , 

3i 

rr 

ax2 

b ^ 

3f 

3i 

0 

a 2, 

3X 

axaf 

3i 

2 

a2, 


E (X )] [ f - E (f)]) 
3i 32 | 


- E (f)]2| 


9f aX3f 


E |[ X - E (X)] [ f - E (f )]2| 


M3<'» = 


{^f „3(X. * ^3(f) 

it X - E ( X 11 ^ [ t - E ( t )lj 
^ E |[ X - E (X)l [ r - E 


3 1 \ 3 



= P 4 (') 


3, \4 /S' \ ^ 

^4 (“afl ^4 



ax 

ai\3 

ar 


E |[X- E(X)]2 [ f- e (f)]| 
+ 4 — ^ E |l X- E(X)] [ f- E (r)]^} 

ax \af/ ' ' 
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The following equations define the various terms in the inclination 
statistics relationships shown on the previous pages- 


E(x) = expected value of x = mean value of x 
AX 

= 

n 


standard deviation of x 


AX- 


1=1 


n-1 


- E(X. 


e|[X — E(X)]*^[Y — E{Y)]^| = expected value of 

[X, - E(X)] '^[Y, - E(Y)]^ 

= _i=2 

n 


/i|^(X) = kth moment of X 
E" [AX, - E(X)] 

= _ i =_1 

n 


Differentiating the relationship i = cos“l (cos X sin f), the following 
partial derivatives are evaluated at the nominal conditions for the desired 
inclination statistics. 


3i sin X sm f 

dX sm 1 


32 , 


cos X sin f sm i — 


in^ f sm^ X cos X 


sm I 


8 - 


3X2 


sm 



di 


- cos X cos f 


sin I 


32 , 


cos X sin f sin i — 


3 9 

COS X cos f sm f 


sin I 


sm 


9^1 

9X3f 


sm X COS f sm I + 


2 2 

sm X sin cos X cos j 
sm I 


sm^ I 


3*5 User Instructions 


MIED utilizes both fixed field and a modified FORTRAN NAMELIST for data 
input. The flight experience data base is input via fixed field format since 
the data base normally does not change. The remaining data, which normally 
changes with each program execution, is input via NAMELIST for flexibility of 
input format. 

Data units are feet, degrees and seconds unless otherwise noted. 

Subsequent data cases are allowed by providing NAMELIST data only. All 
input data is retained on subsequent cases. 

Sample data cases are provided in Appendix A. 

3.3*1 Fixed Field Format - The flight data base is input via fixed field 
format. For each flight the following parameters are input at orbit 
insertion. 


o flight identification 

o option specifying usage of latitude, longitude and azimuth errors. 
0 nominal altitude 
0 nominal latitude 
0 nominal flight azimuth 

o observed errors in altitude, velocity, flight path angle, flight 
azimuth, latitude, longitude. 

This flight data base, consisting of the parameters shown above, was developed 
specifically for orbit accuracy calculations for the Scout launch vehicle. 
Although the use of MIED is not limited to the Scout vehicle, MIED was 
designed to accept this data base without change even though only nominal 
latitude, azimuth and observed errors in latitude, longitude and azimuth are 
used in the MIED calculations. If it becomes desirable to create a new data 
base, the unused parameters need not be included; but the parameters which are 

used must remain within the same column ranges as shown on the next page. 

Fixed field data must either be right justified in the field or 
contain a decimal point. 

The flight experience data base is input in the order shown on the 
next page for each flight. Data may be input for up to 100 flights. 

Preceding the flight data, one title card is used for identification purposes 

only. Following the flight data, an END card beginning in column 1 is used to 
specify the end of the flight data. 
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Fixed Field Input Definition 


Column Range Definition 

1-5 Flight Number 

6-9 Flag 

= 1 if flight IS to be used for latitude and 
longitude statistics only. 

= 2 if flight IS to be used for latitude, longitude 
and azimuth statistics. 

= 3 if* flight IS to be used for azimuth statistics 
only. 

= 4 if flight IS not to be used for statistics. 
10-16 Nominal altitude, n.mi. 

17-23 Altitude error, n.mi. 

24-30 Velocity error 

31-37 Path angle error 

38-44 Azimuth error 

45-52 Latitude error 

53-60 Longitude error 

61-70 Nominal latitude 

71-80 Nominal azimuth 


3.3»2 NAMELIST Format - A modified FORTRAN NAMELIST is used for inputting 
data to MIED. NAMELIST is used because of its readability and simplicity of 
inputting data. The following rules apply to NAMELIST input to MIED: 

!• First card of a data group or case is $INPUTD beginning in 
column 2. Blanks are not allowed. 

2. Last card of a data group or case is $END beginning in column 2. 
Blanks are not allowed. 

3- Blanks may not be used within names but may be used elsewhere. 

4* Variable names are followed by an equal sign which is followed 

by a value which is followed by a comma, e.g., AZIMUTH = 181.5 > 

5« Only columns 2-72, inclusive, are used. 
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6. Titling information may be input by the appropriate title names, 
e.g., TITLEl- THIS CASE PROVIDES INCLINATION ERROR STATISTICS 

7* Any number of names and values may be on a single card or line* 

Definitions of specific NAMELIST inputs to MIED are shown below. Default 
values are shown when they are set by the program prior to reading input 
data. 


NAMELIST Input Definitions 


AZMUTH Nominal geocentric inertial flight azimuth. 

(O. built-in) 

ICOR Non-zero value results in calculation and 

and inclusion of correlation terms in the 
inclination error statistics. (l. built-in) 

IMEAN Non-zero value results in calculation and 

inclusion of mean values in latitude, longitude 
and azimuth in the inclination error 
statistics. (l. built-in) 

LAT Nominal geocentric latitude. (O. built-in) 

SIGAZ If zero, the standard deviation of azimuth is 

calculated from the flight experience sample 
set (O. built-in). If non-zero, this value is 
used for the standard deviation in azimuth. 

TITLEl Title information printed at top of each page of 

TITLE2 output. 72 characters maximum. 


3 .4 Output Description 


Both the fixed field and NAJ-IELIST input data are listed verbatim on the 
output listing (shown in Appendix A on pages A1 and A2) . These lists provide 
a quick check of the input data for format correctness and validity. 

Following the input data lists, output pages 1 and 2 contain information 
related to the calculation of latitude error statistics. Page 1 shows those 
flights with known latitude and longitude errors and the resulting errors in 
range and crossrange. These range and crossrange errors are along and 
perpendicular to the nominal azimuth of each flight. 
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Output page 2 contains adjusted errors in latitude, longitude and 
inclination for each sample on page 1. The range - crossrange errors of each 
flight on page 1 are converted to latitude - longitude errors for the nominal 
flight azimuth input by NAMELIST and shown on page 2. These adjusted errors 
provide a consistent set of sample errors from which to calculate error 
statistics involving latitude and longitude. The azimuth errors shown on page 
2 are those observed on the flights as input. 

Output page 5 contains flight observed azimuth errors. This sample set is 
used to calculate error statistics involving azimuth. 

Output page 4 contains self-explanatory parameters of the inclination 
error statistics calculations. 

Sample data cases are presented in Appendix A. 



TABLE 3.1 


VEHICLE 


S-I3fr 

S-131 

S-133 

S-139 

S-i^O 

S-IA? 

S-1A3 

S-1A3 

S-l^t 

S-1A7 

S-1^0 

S-1^^9 

S-150 

S-l^A 

S-153 

S-155 

S-196 

S-l^7 

S-158 

S-16? 

S-IM 

S-165 

S-167 

S-172 

S-169 

S-17f^ 

S-17^ 

S-179 

S-173 

S-177 

S-lEO 

S-163 

S-1E3 

S-lB<r 

S-1P2 

S-17C 

s-ieB 

S-181 

S-178 

S-190 

s-iee 

S-191 

S-1B6 

S-18Q 

S-l‘»7 

$-194 

S-1Q5 


SOTEi 


SCOUT FLIGHT EXPERIENCE 
AT ORBIT INSERTION 


ERROR FRON PREDICTED 


NSERTICN 

altitude 

ALTITUDE 

VELOCITY 

N,MI, 

N.MI. 

FPS 

56? .e 

9.7 

-58.5 

M2 .5 

6.9 

-280.0 

383 .b 

3. 2 

-107.9 

404 .3 

9.9 

-29.2 

506 ,8 

-.2 

67.2 

4B7.1 

4.1 

142.9 

4fl? .2 

4.2 

112.3 

200.1 

-3.5 

15.0 

495 .2 

-3.7 

-24.3 

350.6 

3,8 

11.1 

197. 7 

2.8 

-31.4 

576.3 

0.0 

-38.6 

178.5 

1.0 

-31.1 

577 .0 

-. 7 

-57.5 

117. 1 

-.5 

-56.9 

279.4 

-5.5 

94.3 

5*^?. 7 

2.0 

-53.4 

565 .9 

2.8 

-34.4 

235.3 

-2.1 

8.0 

531.2 

-1.7 

-29.1 

190.5 

-5.8 

17.0 

376.8 

-2.3 

72.4 

147.2 

.4 

30.7 

21*. 0 

-5.2 

-91.1 

214.7 

-1.3 

-35.0 

588.3 

-10.1 

2.0 

169.9 

-3.1 

-38.2 

294,4 

-.7 

7.6 

115.5 

4.5 

-93.0 

323.3 

.6 

-79.4 

486.0 

4.7 

-202.0 

12J .0 

-.4 

-129.0 

297 .0 

-9.2 

47.0 

263.5 

6.7 

-8.0 

449.6 

8.C 

-177.0 

299.7 

0.0 

-29.5 

151.2 

-3.5 

70.5 

129.5 

-3.4 

84.2 

604.4 

14.5 

-220.2 

123.1 

2.4 

58.2 

405.0 

-2.9 

149.5 

215 .2 

• 6 

45.7 

124 .2 

-4.2 

3.9 

213 .5 

-5.0 

61.1 

270.2 

4.0 

39.1 

271 .0 

.6 

14.1 

194.8 

3.4 

-25.8 


PATH 


ANGLE 

AZINUTH 

LATITUDE 

OEG 

DEG 

DEG 

.77 

• 28 


-.36 

-.15 


.15 

.77 


• 05 

.19 


-.37 

-.91 


-.86 

-.29 


.14 

-.28 


-.75 

.48 


-.63 

0.00 


-.98 

-.24 


-.17 

-.55 


• 16 

-1.18 


-.63 

-.14 

• 0638 

.12 

.26 

-.0473 

-.63 

-.13 


-.06 

.18 

-.1042 

-.05 

-.43 

-.0108 

-.15 

-.73 

-.0429 

-.17 

.70 

-.0729 

-.41 

-.07 

-.2190 

-.46 

-1.08 

-.0774 

-.72 

-1.36 

-.0345 

-.57 

-.25 

-.0300 

-.13 

-.93 

-.0801 

-.13 

• 31 

-.0152 

1.01 

• 02 

-.0233 

-.17 

-.69 

-.0656 

-.07 

.11 

-.213 5 

-.08 

-1.05 


-.63 

-.54 

• 2250 

-.24 

.06 

• 2448 

.49 

-.12 

• 1882 

-;47 

0.00 

• 0015 

• 26 

• 32 

-.1355 

-.02 

.13 

.2266 

-.76 

-.60 

.2533 

1.14 

.34 

-.0743 

-.41 

-.28 

-.0777 

-.03 

.17 

.3180 

.03 

.07 

-.0017 

-.47 

-.62 

-.1845 

-.56 

-.20 

-.0857 

-.26 

.69 

-.0147 

-.05 

• 28 

• 0101 

.15 

0.00 

• 0601 

.03 

-.43 

• 1489 

-.09 

.77 

• 0623 


LONGITUDE 

DEG 


• 0339 
-.0361 

• 0115 
.3134 

.3543 ♦ 

-.0749 

• 4662 
.0137 

• 1321 
.0690 
.0344 
.0477 
.1267 

-.0159 

• 0215 

.0366 

-.1021 

• 7061 
.0170 

-.0692 

.1261 

-.0426 

-.0523 

• 0662 

• 055C 
.0943 
.1735 

-.2036 

-.5603 

-.3837 

• 1696 

• 0124 
-.1025 


LATITUDE, LONGITUDE ERRORS NOT SHOWN ARE NOT AVAILABLE 
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4 . 0 SUBROUTINE DESCRIPTIONS 


This section provides a brief explanation of each subroutine of MIED, 


4*1 MIED (Main Program) 


Data IS initialized - sample error sets are placed in appropriate arrays - 
partial derivatives are evaluated at nominal conditions - ^range and 
crossrange errors are calculated and converted to latitude and longitude 
errors - inclination error statistics are calculated* 


4*2 CARDS 


Subroutine CARDS reads fixed field data in alphanumeric format, writes the 
data on the output file as read, and writes the data on unit 8 for subsequent 
reading by the main program in floating point format. CARDS also counts the 
number of samples input and writes error messages if there are no samples or 
if they exceed the maximum of 100. 


4-5 CORSIG 


Subroutine CORSIG calculates statistical means, standard deviations and 
correlation coefficients from a sample error set of two variables. 


4.4 INPUT 


Subroutine INPUT reads input data cards in a modified NAMELIST format. 
Titling information on title cards are placed in appropriate arrays. 
Non-title cards are written on unit 8 for a NAMELIST read from the main 
program. 


4*5 STAT 


Subroutine STAT calculates statistical expected values and moments from a 
sample error set of two variables. 
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5.0 PROGRAM CODING 


This section presents details about the program coding. Included are 
flowcharts of each subroutine, FORTRAN listings of each subroutine and 
definitions of the FORTRAN variables. The information presented in this 
section is intended to be helpful in developing a thorough understanding of 
MIED and in making modifications to the program. 


5 . 1 Subroutine Flowcharts 


Flowcharts are presented in Figures 5*1 through 5«5* The flowchart 
conventions used are defined in Section 2.0 of this report. 


5.2 FORTRAN Listings 


MIED IS coded in FORTRAN IV, Reference (2), on the CDC CYBER 175 computer 
with the NOS/BE 1.4 operating system. Listings of the FORTRAN coding are 
presented in Appendix B. 


5.3 FORTRAN Variable Definition 


Definitions of the FORTRAN variables are presented below. This 
information is normally used only when making modifications to the program. 


Definition of FORTRAN Variables 
Main Program 


Variables 

ALPHA3 

ALPHA4 

AZMUTH 

AZR 

CAZ 

CEAST 

CLAT 


Definition 

Skewness of inclination 
Kurtosis of inclination 
Input azimuth, deg 
AZMUTH, rad 
Cosine of AZMUTH 

Easterly component of position error, ft 
Cosine of LAT 
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CNORTH 

COV 


CZ 

DAZ 

DOR 

DEG 

DI 

DLAT 

DLG 

DLONG 

DLT 

DR 

DZ 

El 

ELZ, ELZ2, ELZ3 
EL2Z, EL2Z2, EL5Z 

FTNM 

ICODE 

I COR 

IMEAN 

INCL 

INOLR 

IPAGE 

IVEH 

IVEH2 

IWORD 

KASE 

LAT 


Northerly component of position error, ft 

Covariance matrix of latitude, longitude, 
azimuth, deg 

Cosine of azimuth of flight sample 

Azimuth error of input sample data, deg 

Crossrange error, nm 

Degrees per radian 

Inclination error, deg 

Adjusted latitude error, deg 

Longitude error of input sample data, deg 

Adjusted longitude error, deg 

Latitude error of input sample data, deg 

Range error, nm 

Azimuth error of input sample data, deg 

Expected value of inclination, deg 

Expected values of correlation of latitude 
and azimuth 

Feet per nautical mile 
Option of input sample data 
Input correlation option 
Input mean value option 
Inclination, deg 
Inclination, rad 
Page number 

Identifier of input sample data 

Identifier of input sample data. 

Applies to latitude, longitude error 
samples only 

Word 1 of input sample data 

Case number 

Input latitude, deg 
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LATR 

LGMEAN 

LT 

LMEAN 

MEAN I 

M3I 

M3L 

M3Z 

M4I 

M4L 

M4Z 

NFLT 

NFLTX 

NFLT2 

NI4AX 

NSAMP 

PIL 

PILZ 

PIL2 

PIZ 

PIZ2 

RE 

RHOLL 

RHOLOZ 

RHOLZ 


LAT , rad 

Mean value of longitude, deg 

Latitude of input sample data, deg 

Mean value of latitude, deg 

Mean value of inclination, deg 

Third moment of inclination error, deg3 

Third moment of latitude error, deg3 

Third moment of azimuth error, deg3 

Fourth moment of inclination error, deg4 

Fourth moment of latitude error, deg4 

Fourth moment of azimuth error, deg^ 

Number of samples 

Number of samples or zero if 
SIGAZ IS input* non-zero 

Number of samples. Applies to 
latitude, longitude error samples only 

Dimension of error samples 

Number of input sample cards read 

Partial derivative of inclination to latitude 

Partial derivative of inclination to latitude 
and azimuth 

Second partial derivative of inclination to 
latitude 

Partial derivative of inclination to azimuth 

Second partial derivative of inclination to 
azimuth 

Earth radius, ft 

Correlation coefficient of latitude to 
longitude 

Correlation coefficient of longitude to 
azimuth 

Correlation coefficient of latitude to 
azimuth 
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SAZ 

SI 

SIGAZ 

SIGI 

SIGL 

SIGLON 

SIGZ 

SIGZET 

SI2 

SLAT 

SZ 

VARI 

WORD 

ZETA 

ZM 

ZMEAN 


Sxne of AZMUTH 
Sine of inclination 

Input value of standard deviation of azimuth, 
deg 

Standard deviation of inclination, deg 

Standard deviation of latitude, deg 

Standard deviation of longitude, deg 

Standard deviation of azimuth, deg 

Standard deviation of azimuth of latitude, 
longitude error samples only, deg 

Square of sine of inclination 

Sine of LAT 

Sine of azimuth of an individual sample 

Variance of inclination, deg2 

Array of sample error set data 

Array of azimuths of the sample set 

Mean value of azimuth of the latitude, 
longitude error samples only 

Mean value of azimuth of all error samples 
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Figure 5.1 

FLOWCHART OF MAIN PROGRAM MI ED 



CALCULATE EXPECTED 
VALUES AND STATISTICAL 
MOMENTS OF LATITUDE 
AND AZIMUTH 


CALCULATE EXPECTED 
VALUE OF INCLINATION 


CALCULATE VARIANCE, STANDARD 
DEVIATION, THIRD AND FOURTH 
MOMENTS, SKEWNESS AND KURTOSIS 
OF INCLINATION 


TRINT RESULTS 
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Figure 5.2 

FLOWCHART OF SUBROUTINE CARDS 



20 - 





Figure 5.3 

FLOWCHART OF SUBROUTINE INPUT 







Figure 5.4 


FLOWCHART OF SUBROUTINE CORSIG 
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Figure 5.5 

FLOWCHART OF SUBROUTINE STAT 



CALCULATE 
MEAN VALUES 



ACCUMULATE 
SUMMATION OF 
ERRORS TO 3RD 
& 4TH POWER 


SUBTRACT MEAN 
VALUES FROM 
SAMPLE ERRORS 


DO FOR ALL 
SAMPLES 


END LOOP 


IF 

ICOR = 0 


DIVIDE SUMMATIONS 
BY NUMBER OF 
SAMPLES 


ACCUMULATE 
SUMMATION OF 
EXPECTED VALUES 


RETURN 
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6.0 OTHER APPLICATIONS 


Various applications can be made of the information provided in this 
report other than that which computer program MIED was designed to do. As 
described in Section 5*2, the orbital inclination relationship was applied to 
the general expressions for moments of a system of components of correlated 
variables. Since these expressions are general, other relationships of 
correlated variables can be applied in like manner to yield the expected 
value, variance, standard deviation and the third and fourth moments. 

The concept of using flight experience errors to predict errors of future 
flights can also be applied to flight parameters other than those presented in 
this report. An example of another application is Reference (2), which uses 
flight experience to predict errors in the trajectory and orbital parameters. 

Additionally, subroutines STAT, CORSIG and INPUT are written in general 
form and can be used in other computer programs. 
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APPENDIX A 


Sample Data Cases 


45 

FLIGHT EXPEDIENCE 






S-115 

3. 

395.3 

18.7 

-27.8 0.24 

0.58 




S-118 

3. 

396.4 

0.0 

-129.3 -0.03 

0.78 




S-120 

3. 

397.5 

8.5 

-115.6 0.14 

-0.003 




S-113 

3. 

223.1 

9.7 

2.5 0.74 

-0.330 




S-122 

3. 

322.1 

2.4 

-357.7 0.25 

0.114 




S-127 

3. 

152.0 

9.7 

-108.0 0.42 

-0.44 




S-125 

3. 

479.8 

-1.2 

-122.7 0.33 

0.455 




S-134 

3. 

538.2 

6.1 

-114.5 -0.38 

-0.190 




S-12 3 

3. 

495.2 

-8.7 

146.0 0.39 

-0.317 




S-133 

3. 

2i^7.7 

9.5 

-62.4 0.49 

-0.061 




S-135 

3. 

285.0 

4.8 

36.4 -0.17 

-0.382 




S-137 

3. 

115.4 

0.5 

121.4 -0.41 

0.398 




S-136 

3. 

562.8 

9.7 

-58.5 0.766 

0.277 




S-131 

3. 

612.5 

6.9 

-280.0 -0.362 

-0.151 




S-138 

3. 

383.8 

3.2 

-107.9 0.153 

0.77 




S-139 

3. 

404.3 

8.9 

-29.2 0.050 

0.188 




S-140 

3* 

506.8 

-0.2 

67.2 -0.372 

-0.906 




S-142 

3« 

487.1 

4.1 

142.9 -0.063 

-0.29 




S-143 

3. 

482.2 

4.2 

112.3 0.135 

-0.28 




S-145 

3. 

200.1 

-3.5 

15.0 -0.752 

0.475 




S-146 

3. 

495.2 

-3.7 

-24.3 -0.633 

-0.004 




S-147 

3. 

350.6 

3.8 

11.1 -0.975 

-0.244 




S-148 

3. 

197.7 

2.8 

-31.4 -0.168 

-0.545 




S-149 

3« 

576.3 

0.0 

-38.6 0.161 

-1.180 




S-150 

2. 

178.5 

1.0 

-31.1 -0.631 

-0.136 0.0638 

0.0339 

20.7161 

171.556 

S-154 

2. 

577.0 

-0.7 

-57.5 0.124 

0.260 -0.0473 

-0.0361 

12.1707 

100.000 

S-153 

3« 

117.1 

-0.5 

-56.9 -0.632 

-0.134 




S-155 

2. 

279.4 

-5.5 

94.3 -0.055 

0.180 -0.1042 

0.0115 

15.0825 

169.640 

S-156 

2. 

582.7 

2.0 

-53.4 -0.049 

-0.430 -0.0108 

0.3134 

12.0070 

180.000 

S-157 

2. 

565.9 

2.8 

-34.4 -0.151 

-0.727 -0.0429 

0.3543 

12.1900 

180.001 

S-158 

2. 

235.3 

-2.1 

8.0 -0.166 

0.701 -0.0729 

-0.0749 

18.2632 

179.999 

S-162 

2* 

561.2 

-1.7 

-29.1 -0.405 

-0.067 -0.2190 

0.4662 

12.2504 

160.000 

S-161 

2* 

190.5 

-5.0 

17.0 -0.459 

-1.083 -0.0774 

0.0137 

22.0366 

186.853 

S-165 

2. 

376.8 

-2.3 

72.4 -0.716 

-1.361 -0.0345 

0.1321 

13.4573 

171.761 

S-167 

2. 

147.2 

0.4 

30.7 -0.569 

-0.250 -0.0300 

0.0690 

23.1376 

184.350 

S-172 

2. 

216.0 

-5.2 

-91.1 -0.130 

-0.929 -0.0801 

0.0344 

20.8663 

175.719 

S-16 9 

2. 

214.7 

-1.3 

-35.0 -0.13 

0.31 -0.0152 

0.0477 

17.7236 

193.313 

S-176 

2. 

568.3 

-10.1 

2.0 -l.OI 

0.02 -0.0233 

0.1267 

11.9576 

180.000 

S-174 

2. 

169.9 

-3.1 

-38.2 -0.164 

-0.608 -0.0656 

-0.0159 

35.5464 

103.438 

S-175 

2. 

294.4 

-0.7 

7.6 -0.072 

0.109 -0.2135 

0.0215 

-2.7201 

88.955 

S-173 

3. 

115.5 

4.6 

-92.4 -0.074 

-1.043 




S-177 

2. 

323.3 

0.6 

-79.0 -0.63 

-0.54 0.2250 

0.0386 

23.3147 

137.246 

S-100 

2. 

466.0 

4.7 

-202. -0.24 

0.06 0.2448 

-0.1021 

23.5077 

135.496 

S-163 

2. 

120.0 

-0.4 

-129. 0.49 

-0.12 0.1682 

0.7061 

-3.3437 

90.865 

S-18 3 

2. 

297.0 

-9.2 

47.0 -0.47 

0. 0.0015 

0.0170 

15.4234 

172.737 


END 


A1 



lINPUTO 

TITLEl-TEST CASE 

TITLE2-WITH MEAN WITH CORRELATION 

IMEAN«1# ICOR-1* 
LAT«15.263»AZMUTH»187.27A» 

SEND 


DATE IS 12/15/80 
TIME IS 12. 36. AS 
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VQUGHT CORPORATION 
PROGRAN NIEO 
CASE 1 PAGE 1 

TEST CASE 

WITH KEAN WITH CORRELATION 


FLIGHT EXPERIENCE ERRORS 



ONLY 

SAMPLES WITH 

LATITUDE, 

LONGITUDE 

ERRORS 



VEHICLE 

LATITUDE 

LONGITUDE 

AZIMUTH 

RANGE 

C-RANGE 



DEG 

DEG 

DEG 

NMI 

NMI 

1 

S-150 

• 06A 

.034 

-.136 

-3.513 

2.448 

2 

S-154 

-.047 

-.036 

.260 

2.843 

-2.121 

3 

S-155 

-.104 

• 012 

.180 

6.261 

-.470 

4 

S-156 

-.Oil 

• 313 

-.430 

.649 

18.420 

5 

S-157 

-.043 

.354 

-.727 

2.578 

20.816 

6 

S-158 

-.073 

-.075 

.701 

4.362 

-4.275 

7 

S-162 

-.219 

.466 

-.067 

13.164 

27.364 

e 

S-161 

-.077 

.014 

-1.083 

4.479 

1.470 

9 

S-165 

-.035 

.132 

-1.361 

3.159 

7.345 

10 

S-167 

-.030 

• 069 

-.250 

1.509 

3.940 

11 

S-172 

-.080 

.034 

-.929 

4.945 

1.567 

12 

S-169 

-.015 

• 048 

• 310 

.260 

2.866 

13 

S-176 

-.023 

.127 

.020 

1.401 

7.450 

l<r 

S-17A 

-.066 

-.016 

-.688 

.160 

-4.016 

15 

S-175 

-.214 

.022 

.109 

1.057 

-12.854 

16 

S-177 

.225 

.039 

-.540 

-8.484 

10.746 

17 

S-180 

.24 5 

-.102 

.060 

-14.439 

6.301 

16 

S-163 

.188 

.706 

-.120 

42.194 

11.951 

19 

S-183 

.002 

.017 

0.000 

• 035 

• 989 
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V0U6HT corporation 
PROGRAM MIEO 
CASE 1 PAGE 2 

TEST CASE 

WITH MEAN WITH CORRELATION 


ADJUSTED ERRORS 

ONLY SAMPLES WITH LATITUDE# LONGITUDE ERRORS 
AZIMUTH- 187,274 DEG 



LATITUDE 

LONGITUDE 

AZIMUTH 

INCLINAT 


OEG 

DEG 

DEG 

OEG 

1 

.053 

• 050 

-.136 

-.133 

2 

-.042 

-.042 

.260 

.252 

3 

-.103 

-.022 

.180 

.177 


-.050 

.314 

-.430 

-.413 

5 

-.086 

.350 

-.727 

-.698 

6 

-.063 

-.083 

.701 

.678 

7 

-.275 

.440 

-.067 

-.056 

8 

-.077 

.015 

-1.083 

-1.042 

Q 

-.068 

.119 

-1.361 

-1.311 

10 

-.033 

• 064 

-.250 

-.240 

11 

-.085 

.C16 

-.929 

-.893 

12 

-.010 

• 048 

.310 

.299 

13 

-.039 

• 124 

• 020 

• 021 

14 

• 006 

-•069 

— • 688 

-.664 

15 

.010 

-.222 

.109 

.105 

16 

• 117 

• 202 

-.540 

-.524 

17 

.225 

.139 

.060 

• 050 

18 

-.722 

.112 

-.120 

-.092 

19 

-.003 

• 017 

0.000 

• 000 


COVARIANCE 

MATRIX based ON 19 

FLIGHTS 

LATITUDE 

-LONGITUDE-AZIMUTH 

-OEG 

•035163 

-.005174 

•000343 

-•005174 

.025925 

-.020693 

•000343 

-.020693 

•281815 
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VQUGHT CORPORATION 
PR06RAN NIEO 
CASE 1 PAGE 3 

TEST CASE 

WITH HEAN WITH CORRELATION 


A7IHUTH DEVIATIONS(OEG) 
FLIGHT EXPERIENCE 


1 

S-115 

.580 

2 

S-118 

.700 

3 

S-120 

-.003 

‘f 

S-113 

-.338 

5 

S-122 

.114 

6 

S-127 

-.440 

7 

S-125 

.455 

8 

S-134 

-.190 

9 

S-123 

-.317 

10 

S-133 

-.061 

11 

S-135 

-.382 

12 

S-137 

.398 

13 

S-136 

.277 

1^ 

S-131 

-.151 

15 

S-138 

.770 

16 

S-139 

.180 

17 

S-140 

-.906 

10 

S-142 

-.290 

19 

S-143 

-.280 

20 

S-145 

.475 

21 

S-146 

-.004 

22 

S-147 

-.244 

23 

S-140 

-.545 

24 

S-149 

-1.180 

25 

S-150 

-.136 

26 

S-154 

• 260 

27 

S-153 

-.134 

28 

S-155 

.180 

29 

S-156 

-.430 

30 

S-157 

-.727 

31 

S-158 

.701 

32 

S-162 

-.067 

33 

S-161 

-1.083 

34 

S-165 

-1.361 

35 

S-167 

-.250 

36 

S-172 

-.929 

37 

S-169 

.310 

38 

S-176 

.020 

39 

S-174 

-.608 

40 

S-175 

.109 

41 

S-173 

-1.043 

42 

S-177 

-.540 

43 

S-180 

• 060 

44 

S-163 

-.120 

45 

S-183 

0.000 
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VOUGHT CORPORATION 
PROGRAM HIED 
CASE 1 PAGE A 

TEST CASE 


WITH MEAN WITH CORRELATION 

NOHINAL VALUES 

inclination# deg 

97.0161 

LATITUDE# deg 

15.2630 

AZIHUTH# deg 

187.274 

INCLINATION STATISTICS 

MEAN# DEG 

-.151370 

VARIANCE# DEG 

.248552 

STANDARD DEVIATION, DEG 

.498550 

THIRD MOMENT 

-.292903E-06 

SKEWNESS 

-.444597 

FOURTH MOMENT 

.142631E-07 

KURTOSIS 

2.48812 

MEAN VALUES 

LATITUDE* DEG 

-.655065E-01 

LONGITUDE* DEG 

.827947E-C1 

AZIHUTHC19 FLIGHTS)* DEG 

-.246895 

AZINUTH(A5 FLIGHTS)* DEG 

-.159156 

STANDARD DEVIATIONS 

LATITUDE* DEG 

. 187518 

LONGITUDE* DEG 

.161012 

AZINUTH(19 FLIGHTS)* DEG 

.530862 

AZIMUTH(A5 FLIGHTS)* DEG 

.517076 

CORRELATION COEFFICIENTS 

LATITUDE-AZIMUTH 

•344954E-02 

LATITUDE-LONGITUDE 

-.171356 

LONGITUDE-AZIMUTH 

-.242089 


TERMS OF THE INCLINATION 

STATISTICAL EQUATIONS 


( 1-INCLINATION} L-LATITUDE) Z 

•AZIMUTH) 

PARTIAL DERIVATIVES 

EXPECTED VALUES 


UNITLESS 


RADIANS 

PIL 

-.335827E-01 

ELZ 

.104602E-06 

PIL2 

-.122931E+00 

EL2Z 

•232265E-07 

PIZ 

.96A183E+00 

ELZ2 

.141576E-07 

PIZ2 

-.865801E-02 

EL2Z2 

•112616E-09 

RILZ 

-.267087E+00 

EL3Z 

-.153411E-09 



ELZ3 

.144715E-09 


MOMENTS 

VARIANCE TERMS 


RADIANS 


RADIANS 

M3L 

-.721875E-07 

TERMl 

.120801E-07 

M4L 

•955A17E-09 

TERM2 

.757152E-04 

M3Z 

-.325381E-06 

TERM3 

-.677401E-08 

M4Z 

.165228E-07 

TERM4 

-.298017E-09 



TERMS 

.271625E-08 



TERM6 

.411645E-11 



TERM7 

-.275300E-08 



TERM8 

.416660E-09 



TERM9 

-.729172E-08 
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WITH CORRELATION 


SINPUTD 

TITLE2- NO MEAN 
IMEAN-0* ICOR-l» 
SEND 


DATE IS 12/15/80 
time is 12.36.A7 
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TEST CASE 

NO MEAN WITH CORRELATION 


VOUGHT CORPORATION 
PROGRAM MIED 
CASE 2 PAGE 1 


FLIGHT EXPERIENCE ERRORS 



ONLY 

SANPLES WITH 

LATITUDE# 

LONGITUDE 

ERRORS 



VEHICLE 

LATITUDE 

LONGITUDE 

AZIMUTH 

RANGE 

c-ran5e 



DEG 

DEG 

DEG 

NMI 

NHI 

1 

S-150 

• 064 

.034 

-.136 

-3.513 

2.448 

2 

$-154 

-.047 

-.036 

.260 

2.843 

-2.121 

3 

S-155 

-.104 

.012 

• 180 

6.281 

-.470 


$-156 

-.011 

.313 

-.430 

.649 

18.420 

5 

$-157 

-.043 

• 354 

-.727 

2.578 

20.816 

6 

$-158 

-.073 

-.075 

.701 

4.382 

-4.275 

7 

$-162 

-.219 

• 466 

-.067 

13.164 

27.384 

8 

$-161 

-.077 

• 014 

-1.083 

4.479 

1.470 

9 

$-165 

-.035 

.132 

-1.361 

3.159 

7.345 

10 

$-167 

-.030 

.069 

-.250 

1.509 

3.940 

11 

$-172 

-.080 

.034 

-.929 

4.945 

1.567 

12 

$-169 

-.015 

.048 

.310 

• 260 

2.868 

13 

$-176 

-.023 

.127 

.020 

1.401 

7.450 

14 

$-174 

-.066 

-.016 

-.688 

.160 

-4.016 

15 

$-175 

-.214 

.022 

.109 

1.057 

-12.854 

16 

$-177 

.225 

.039 

-.540 

-8.484 

10.746 

17 

$-180 

.245 

-.102 

.060 

-14.439 

6.301 

18 

$-163 

• 186 

.706 

-.120 

42.194 

11.951 

19 

$-183 

.002 

.017 

0.000 

• 035 

.989 
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TEST CASE 

NO hean with correlation 


V0U6HT CORPORATION 
PROGRAM HIED 
CASE 2 PAGE 2 


ADJUSTED ERRORS 

ONLY SAMPLES WITH LATITUDE* LONGITUDE ERRORS 
AZIMUTH* 187. 27A DEG 



LATITUDE 

LONGITUDE 

AZIMUTH 

INCLINATION 


DEG 

DEG 

DEG 

DEG 

1 

.053 

.050 

-.136 

-.133 

2 

-.042 

-.042 

• 260 

.252 

3 

-.103 

-.022 

• 180 

.177 

4 

-.050 

.314 

-.430 

-.413 

5 

-.086 

.350 

-.727 • 

-.698 

6 

-.063 

-.083 

.701 

.678 

7 

-.275 

.440 

-.067 

-.056 

8 

-.077 

.015 

-1.083 

-1.042 

Q 

-.068 

.119 

-1.361 

-1.311 

10 

-.033 

.064 

-.250 

-.240 

11 

-.085 

.016 

-.929 

-.893 

12 

-.010 

• 048 

.310 

.299 

13 

-.039 

.124 

.020 

.021 

14 

.006 

-.069 

-.688 

-.664 

15 

.010 

-.222 

.109 

.105 

16 

.117 

.202 

-.540 

-.524 

17 

.225 

.139 

.060 

• 050 

16 

-.722 

.112 

-.120 

-.09 2 

19 

-.003 

.017 

0.000 

• 000 


COVARIANCE MATRIX BASED ON 19 FLIGHTS 
LATITUDE-LONGITUDE-AZIMUTH-OEG 


,039454 

-.010597 

.016517 


-.010597 

.032780 

-.041134 


.016517 

.041134 

.342772 
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TEST CASE 

NO HEAN WITH CORRELATION 


VOUGHT 

CASE 


CORPORATION 
PROGRAM HIED 
2 PAGE 3 


AZIMUTH DEVIATIONS(OEG) 



FLIGHT 

EXPERIENCE 

1 

S-115 

• 580 

2 

S-118 

• 78 0 

3 

S-120 

-•003 

A 

S-113 

-•338 

5 

S-122 

• 114 

6 

S-127 

-.440 

7 

S-125 

• 455 

8 

S-134 

-.190 

9 

S-123 

-•317 

10 

S-133 

-.061 

11 

S-135 

-.382 

12 

S-137 

.398 

13 

S-136 

.277 

14 

S-131 

-.151 

15 

S-138 

.770 

16 

S-139 

.168 

17 

S-140 

-.906 

18 

S-142 

-.290 

19 

S-143 

-.280 

20 

S-145 

.475 

21 

S-146 

-.004 

22 

S-147 

-.244 

23 

S-148 

-.545 

24 

S-149 

-1.180 

25 

S-150 

-.136 

26 

S-154 

• 260 

27 

S-153 

-.134 

28 

S-155 

• 180 

29 

S-156 

-.430 

30 

S-157 

-.727 

31 

S-158 

.701 

32 

S-162 

-.067 

33 

S-161 

-1.083 

34 

S-165 

-1.361 

35 

S-167 

-.250 

36 

S-172 

-.929 

37 

S-169 

.310 

38 

S-176 

.020 

39 

S-174 

-.688 

40 

S-175 

.109 

41 

S-173 

-1.043 

42 

S-177 

-.540 

43 

S-180 

.060 

44 

S-163 

-.120 

45 

S-183 

0.000 
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TFST CASE 

NO MEAN WITH CORRELATION 


VOUGHT CORPORATION 
PROGRAM 1IEO 


CASE 2 PAGE A 


NOMINAL VALUES 

INCLINATION# DEG 
LATITUDE# DEG 

azimuth, deg 

INCLINATION STATISTICS 
MEAN, DEG 
VAPIANCE# DEG 
STANDARD DEVIATION, DEG 
THIRD moment 
SKEWNESS 
FOURTH MOMENT 
KUPTOSIS 
MEAN VALUES 

LATITUDE, DEG 
LONGITUDE, DEG 
AZIMUTHdR FLIGHTS)# DEG 
AZIMUTH(A5 FLIGHTS), DEG 
standard DEVIATIONS 
LATITUDE, DEG 
LONGITUDE, DEG 
AZIMUTHd? FLIGHTS), DEG 
AZIMUTH<A5 FLIGHTS), DEG 
CORRELATION COEFFICIENTS 
LATITUDE-AZIHUTH 
LATITUDE-LQN6ITU0F 
LONGITUDE-AZIMUTH 


97.0161 
15.2630 
187. 27A 

-.1AIA33E-03 

.271297 

.520862 

-.12B516E-05 

-1.7106A 

.269679E-07 

3.9A86A 

0 . 

0 . 

0 . 

0 . 

.198630 

.181052 

.585A67 

.541016 

. 142028 
-.294676 
-.388058 


TERMS OF THE INCLINATION 

STATISTICAL equations 


(I-INCLINATIONj L-LATITUOEJ Z 

•AZIMUTH) 

PARTIAL DERIVATIVES 

EXPECTED VALUES 


UNITLESS 


RADIANS 

PIL 

-.335827E-01 

EtZ 

.503125E-05 

PIL? 

-•122931E+00 

EL2Z 

-.260757E-07 

PIZ 

• 96<tie3E+00 

EIZ2 

-•998382E-07 

PIZ2 

-.865801E-02 

EL2Z2 

.19aS8lE-09 

PILZ 

-.267087E+00 

EL3Z 

.233801E-09 



ELZ3 

• 1M858E-08 


HCIMENTS 

VARIANCE TERMS 


RADIANS 


RADIANS 

M3L 

-.108251E-06 

TERMl 

.135543E-07 

M4L 

•136630E-08 

TERH2 

.828885E-04 

M3Z 


TEPM3 

-.325823E-06 

M4Z 

•314280E-07 

TERM4 

-.446900E-09 



TERNS 

.120553E-07 



TERM6 

-.290289E-10 



TERM7 

.309070E-08 



TERM8 

-.467772E-09 



TERM9 

.514208E-07 
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APPENDIX B 


FORTRAN Code Listings 


PROGRAM NIE0( INPUT# OUTPUT# TAP E 5-IN PUT# TAPEt- OUTPUT# TAPER) 

CALCULATES ORBITAL INCLINATION STATISTICS FROM A SAMPLE SET 
OF AZIMUTH# LATITUDE AND LONGITUDE ERRORS OBTAINED FROM FLIGHT 
EXPERIENCE. EQUATIONS OBTAINED FROM A TAYLOR SERIES EXPANSION, 
REFERENCE# ^STATISTICAL MOOEIS IN ENGINEER INGi«# BY HAHN AND SHAPIRO, 


«#*«**««««««««* FIXED FIELD INPUTS •*♦♦♦♦♦*♦♦** 

WORD 1 • VEHICLE NUMBER (COLUMNS 1-5) 

WORD 2 ■ DATA OPTION (COLUMNS 6-9) 

-1 FLIGHT SAMPLE IS USED FOR LATITUDE, LONGITUDE 
STATISTICS ONLY 

•2 FLIGHT SAMPLE IS USED FOR LATITUDE, LONGITUDE 
AND AZIMUTH STATISTICS 

-3 FLIGHT SAMPLE IS USED FOR AZIMUTH STATISTICS ONLY 
•A FLIGHT SAMPLE IS NOT USED FOR STATISTICS 


WORD 

WORD 

WORD 

WORD 

WORD 

WORD 

WORD 

WORD 

WORD 


3 

A 

5 

6 

7 

8 
9 

10 

11 


NOMINAL ALTITUDE# N.M. 
ALTITUDE ERROR# N.M. 
VELOCITY ERROR# FPS 
PATH ANGLE ERROR# DEG 
AZIMUTH ERROR# DEG 
LATITUDE ERROR# DEG 
LONGITUDE ERROR# DEG 
NOMINAL LATITUDE# DEG 
NOMINAL AZIMUTH, DEG 


(COLUMNS 10-16) 
(COLUMNS 17-23) 
(COLUMNS 2A-30) 
(COLUMNS 31-37) 
(COLUMNS 38-AA) 
(COLUMNS A5-52) 
(COLUMNS 53-60) 
(COLUMNS 61-70) 
(COLUMNS 71-80) 


^tttt********* NAMELIST INPUTS ♦♦♦♦♦♦♦♦♦♦♦♦* 

AZMUTH- NOMINAL GEOCENTRIC INERTIAL AZIMUTH. 

ICOR - NON-ZERO VALUE RESULTS IN INCLUSION OF CORRELATION TERMS IN THE 
INCLINATION STATISTICS. 

IMEAN - NON-ZERO VALUE RESULTS IN CALCULATON OF MEAN VALUES IN 

LATITUOE#LONGITUDE#AZIMUTH. ZERO VALUE ASSUMES ZERO MEANS. 

LAT - NOMINAL GEOCENTRIC LATITUDE. 

SIGAZ - ZERO VALUE RESULTS IN USING INPUT AZIMUTH ERROR SAMPLE SET TO 
CALCULATE STANDARD DEVIATION, NON-ZERO VALUE RESULTS IN USING 
SIGAZ AS STANDARD DEVIATION. 

TITLEl- TITLE PRINTED AT TOP OF EACH PAGE. 72 CHARACTERS MAXIMUM. 
TITLE2- TITLE PRINTED AT TOP OF EACH PAGE. 72 CHARACTERS MAXIMUM. 


IMPLICIT REAL ( A-H, L, M# 0-Z ) 
integer TITLE1#TITLE2 
REAL INCL#INCLR#INCLN 

COMMON /STATS/ EL Z# EL2Z# EL Z2# EL3Z#EL Z3#EL2Z2# M3L# M3Z# MAL# M4Z 
DIMENSION 0R(100)#0CR(100)#DZ(100)#DLAT(10O)#DL0NG(lDO)# 

* OAZ(100)#DLT(100)#DLG(100)#ZETA(100),IVEH(100), IVEH2(100), 

- LT( 100), WORD( 11# 100) ,TITLE1(8 ) #TITLE2(8 )#COV (3#3) #DI (100)# 
A IWORD(IOO) 

DATA TITLE1#TITLE2/16*10H / 

DATA deg/57.29578/# RE / 20925741 . /# FTNM/6076, 1155/ 

NAMELIST /INPUTD/ LAT# AZMUTH# IMEAN# I COR# S IGAZ 
1NCLN(L# A)-ACOS (COS (L/OEG)*SIN(A/DEG) )*0EG 
C 

c 

WRITE(6#70) 

C INITIALIZE 

AZMUTH-0. 
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ICOB-O 

THE«N-0 

l*T»0, 

SIGAZ>0. 

NHAX-100 

KASE-0 

READ FIXED FIELD INPUT DATA 
CALL CARDS(NSANP>NHAX) 

READ! 8« 80) ( IU0RD( J ) » ( WORD ( I> J ) t I >2» 11 ) * J ■!* NS AHP ) 
READ NAMELIST INPUT DATA 
10 CONTINUE 
KASE-KASE+1 

CALL INPUT (TITLE1»TITLE2) 

READ(8>INPUT0) 


PLACE SAMPLE SET DATA IN APPPOPIATE ARRAYS 
NFLT-0 
NFLT2-0 

DD 40 I-1»NSAMP 
IC00E-W0P0(2>II 
GOTO (20>20*30>40)» ICODE 

SELECTS SAMPLES TO BE USED FOR LATITUDEt LONGITUDE STATISTICS 
20 NFLT2-NFLT2*! 

IVEH2(NFLT2)«IWOROII) 

DZ(NFLT2)-W0R0(7>I) 

0LT(NFLT2)-W0RD(8»I) 

0LG(NFLT2)-W0RD(9>1) 

LT (NFLT2)-W0RD(10»I) 

ZETA(NFLT2)-W0R0( 11» I ) 

GOTO (40>30>30«40)* ICODE 
SELECTS SAMPLES TO BE USED FOR AZIMUTH STATISTICS 
30 NFLT-NFLT+1 

IVEH(NFLT)>IWOPO(I) 

DAZ(NFLT)-W0RD(7»I) 

40 CONTINUE 
SET-UP 

AZP-AZHUTH/DEG 

LATR-LAT/OEG 

SAZ>SIN(AZR) 

CAZ-COS(AZR) 

SLAT-SIN(LATP) 

CLAT-COS (LATR) 

INCL-INCLN(LAT*AZMUTH) 

INCLR-INCL/OEG 

SI-SIN(INCLR) 

SI2-SI**2 

NFLTX-NFLT 


COMPUTE PARTIAL DERIVATIVES AT NOMINAL CONDITIONS 
PIL-SLAT*SAZ/SI 

PIL2-ICLAT*SAZ4SI-SAZ*P3*SLAT**2*CLAT/SII /SI2 
PIZ— CLAT4CAZ/SI 

PIZ2-CCLAT4SAZ*SI-CLAT*43*SAZ*CAZ**2/SI)/SI2 
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PILZ-(SLAT*CAZ*SI+SLAT*SAZ**2*CLAT**2*CAZ/SI) /SI2 

C 

C CONVERT LATITUDE-LONGITUDE DEVIATIONS TO P ANGE-CROSSRANGE DEVIATION 

IPAGE-1 

WRITE (6*90) KASE* I PAGE* TIT LEI* TITLE 2 

WRITE(6*100) 

no 50 I-1*NFLT2 

SZ«SIN(ZETA(I)/DEG) 

CZ-COS(ZETA( I) /DEG) 

CN0RTH-RE*DLT(I)/DEG 
CEAST-RE*DLG(I)*COS(LT(I)/DEG)/OEG 
DR (I)»(CNORTHACZ^CEAST*SZ) /FTNM 
DCR(I)-(CN0RTH*SZ-CEAST*CZ)/FTNM 

WRITE (6* 110) I*IVEH2 (I )*DLT(I)*DLG( I)*DZ( I)*DR (I )*DCR (I ) 

50 CONTINUE 

C CONVERT RANGE-CROSSRANGE DEVIATIONS TO LATITUDE-LONGITUDE DEVIATION 

DO 60 I«1*NFLT2 

CNORTH»DR(I)*CAZ+DCR(I)*SAZ 
CEAST«DRII)*SAZ-DCR( I)*CAZ 
DLAT(I)-CNORTH/RE*FTNM*OEG 
nLONG(I)»CEAST/RE*FTNM/CLAT*DEG 
DI(I)-INCLNtLAT*OLAT(I)*AZNUTH*DZtI)) - INCL 
60 CONTINUE 
IPAGE-2 

WRITE (6*90) KASE*IPAGE*TITLE1*TITLE2 
WRITE(6*120) AZMUTH 

WRITE (6* 130) (I*0LAT(I)*DL0NG(I)*DZ(I)*DI(I)*I-ltNFLT2) 


C 


C 

C 


C 


COMPUTE SIGMA LATITUDE* LONGITUDE* AZIMUTH AND CORRELATIONS (FLT. 
PARAMETERS ARE FOR LATITUDE-LONGITUDE SAMPLES ONLY 
CALL CORSIG (0LAT*DZ*NFLT2» SIGL* SIGZET*RHOLZ* 

LTMEAN,ZM*IMEAN*IC0R) 


CALL CORSIG (DLONG* 0Z*NFLT2* S IGLON* S IGZET* RHOLOZ* 
A LGMEAN*ZM*IMEAN*IC0R) 

CALL CORSIG ( DL AT*OLONG* NFLT2* SIGL* SIGLON »RHOLL , 

A LTMEAN*LGMEAN*IHEAN*IC0R) 

COMPUTE COVARIANCE MATRIX 
C0V(1*1)»SIGL**2 
C0V(1*2)»RH0LL*SIGL*SIGL0N 
C0V(1*3)«RH0LZ*SIGL*SIGZET 
COV(2*l)>COV(l*2) 

C0V(2*2)»SIGL0N**2 

COV(2*3)-RHOLOZ*SIGLON«SIGZET 

C0V(3*1)-C0V(1»3) 

C0V(3*2)«C0V(2*3) 

C0V(3*3)»SIGZET*R2 

WRITE (6, 140) NFLT2* ((C0V(1*J)*J»1*3)*I»1»3) 


EXP. ) 


CALCULATE SIGMA AZIMUTH FROM FLIGHT EXPERIENCE( ALL FLIGHTS) 
IPAGE-3 

WRITE (6*90) KASE*IPAGE*TITLE1»TITLE2 
WRITE (6*150) (I*IVEH(I)*0AZ(I)*I»1»NFLT) 

ZMEAN • AZIMUTH MEAN OF ALL FLIGHTS 
CALL CORSIG (OAZ* DAZ*NF LT* S IGZ*DUhl* DUM2* ZHE AN* 0UM3* IME AN. ICOR ) 
IF (SIGAZ.NE.O.) SIGZ-SIGAZ 
IF (SIGAZ.NE.O.) NFLTX-0 
IF (SIGAZ.NE.O.) ZMEAN-0. 
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SI6L-SI6L/0EG 

SIGZ-SIGZ/DEG 


CGHPUTE INCLINATION STATISTICS 

COMPUTE LATITUDE* AZIMUTH STATISTICS 
CALL STAT(DLAT*DZ*NFLT2*IC0R*IMEAN) 

C 

C COMPUTE EXPECTED VALUE DF INCLINATION 

FIl-INCLN(LAT*LTMEAN,AZMUTH+ZMEAN)/nEG 
EI2-0.5PP1L2PSIGLP*2 
EI3.-0.5*PIZ2*SIGZ*^2 
EI4 -pilZ*ELZ 

EI«(EH*EI2*EI3*EIA)*DEG 

MEANI-EI-INCL 

C 

C COMPUTE VARIANCE OF INCLINATION 

TERM1-PIL**2*SIGL**2 
TERM2«PIZ**2PSIGZ**2 
TERM3-2.*PIL*PIZ*ELZ 
TERMA-PILPPIL2PM3L 
TERM5-PIZ*PIZ2*M32 
TERM6»PIL*PIZ2*ELZ2 
TERM7»PIZ*PIL2*EL2Z 
TERM8»2.PPIL*PILZ*EL2Z 
TBRM<)-2.*PIZ*PILZ*ELZ2 

VARI-TERM1<-TEPM2+TERM3+TERM4+TERM5+TERM6*TERM7+TERM8+TERM9 
SI6I«S0RT(VARI J*DE6 
C 

C COMPUTE THIRD MOMENT AND SKEWNESS OF INCLINATION 

M3I»PIL**3*M3L+PIZ**3PM3Z+3.*PIL**2PPIZ*EL2Z+3.*PIL*PIZ**2PELZ2 
ALPHA3-M3I/SQRT( VARI*P3) 
r 

C COMPUTE FOURTH MOMENT AND KURTOSIS OF INCLINATION 

M4I«pILP*4*M4L+PIZ**4*M4Z44.*PILP*3*PIZ*EL3Z+4.*PIL*PIZ**3*ELZ3+6. 
A ♦PIL**2*PIZ**2*EL2Z2 
ALPHA4«M4I/VARIP*2 
C 

VARI«VARI*DEG**2 

SIGL-SIGL*DEG 

SIGZ-SIGZ*OEG 

IPAGE-4 

WRITE (6, 90) KASE,IPAGE*TITLE1,TITLE2 

WRITE 16*160 ) INCL*LAT*AZMUTH*MEANI*VARI*SIGI*M3I*ALPHA3*M4I, ALPHA4 
WRITE (6*170) LTMEAN,LGMEAN,NFLT2*ZM*NFLT,ZMEAN,SIGL»SIGL0N, 

A NFLT2*SIGZET»NFLTX»S1GZ*PHOLZ*RHOLL.RHOLOZ 

WRITE (6* 180) PIL*ELZ*PIL2*EL2Z*PIZ,ELZ2*PIZ2.EL2Z2.PILZ,EL3Z,ELZ3 
WRITE (6* 190) M3L.TERM1.M4L»TEPM2*M3Z*TERM3*M4Z*TFRM4. 

A TERM5*TERM6*TERM7,TERM8»TERM9 

GOTO 10 
C 

70 FORMAT (IHI) 

80 format (A5*F4.0*5F7.0»2F8.0*2F10.0) 

90 FORMAT ( 1H1/T66**V0UGHT CORPORATION*, /T72**PR0GRAM MIEO**/ 

A T66,*CASE *12* PAGE *,11/ 2 ( lOX, 8A10/ ) / ) 

100 FORMAT (26X,*FLIGHT EXPERIENCE ERRORS*,// 
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A 16X»*0NLY samples WITH LATITUDE* LONGITUDE ERRORS**// 

A IIX** VEHICLE LATITUDE LONGITUDE AZIMUTH RANGE** 

A ♦ C-RANGE**/ 

A 21X*3(7X**DEG*)*?(7X**NMI*I/) 

110 FORMAT ( lOX* I2*2X*A5»3X*5F10.3) 

120 FORMAT (31X,*ADJUSTED ERRORS**// 

A 17X**0NLY SAMPLES WITH LATITUDE* LONGITUDE ERRORS**// 

A 29X**AZIMUTH»**F8.3,* DEG*) 

130 FORMAT < A( / )23X**LATITUDE LONGITUDE AZIMUTH INCLINATION**/ 

A 19X*4(7X,*DEG*) //50(18X,I2*AF10.3/)) 

lAO FORMAT (6(/)22X**C0VARIANCE MATRIX BASED ON*, 13** FLIGHTS**/ 

A 25X**LATITUDE-LONGITUOE-AZIMUTH-DEG*, // 


A 

3(15X,3F15.6W ) ) 



150 FORHAT 

( 15X**AZIMUTH OEVIATIONS(OEG)**/ 


A 

20X**FLIGHT EXPERIENCE** 

// 


A 

100(15X*I2»3X* A5*3X*F7.3 

#/)) 


160 FORMAT 

( 10X**N0MINAL VALUES 


♦ » / 

A 

lOX** INCLINATION* DEG 

♦#G15.6/ 

A 

lOX** LATITUDE, DEG 


♦#G15.6/ 

A 

lOX,* AZIMUTH* DEG 


♦»G15*6/ 

A 

10X**INCLINATI0N STATISTICS 

/ 

A 

lOX** MEAN* DEG 


«#G15.6/ 

A 

lOX** VARIANCE* DEG 


*#G15.6/ 

A 

lOX** standard DEVIATION* DEG 

♦»G15*6/ 

A 

lOX** THIRD moment 


♦#G15.6/ 

A 

lOX** SKEWNESS 


♦#G15.6/ 

A 

lOX** FOURTH MOMENT 


♦#615,6/ 

A 

lOX,* KURTOSIS 


♦fG15.6) 

170 FORMAT 

( 10X,*MEAN VALUES 


*9 / 

A 

lOX** LATITUDE* DEG 


♦#G15.6/ 

A 

lOX** LONGITUDE* DEG 


«»G15«6/ 

A 

lOX** AZIMUTH!** 12** 

FLIGHTS)# 

DEG >^#G15.6/ 

A 

lOX,* AZIMUTH!*, 12** 

FLIGHTS)# 

DEG ♦»G15«6/ 

A 

10X**STANDARD DEVIATIONS 

*9 / 

A 

lOX** LATITUDE* DEG 


♦#G15.6/ 

A 

lOX** LONGITUDE* DEG 


♦#G15.6/ 

A 

lOX,* A7IMUTH!**I2,* 

FLIGHTS)* 

DEG t,G15«6/ 

A 

lOX** AZIMUTH!** 12** 

FLIGHTS), 

DEG <»#G15.6/ 

A 

10X**C0RRELATI0N COEFFICIENTS 

♦# / 

A 

lOX** LATITUDE-AZIMUTH 

♦#G15.6/ 

A 

lOX,* LATITUDE-LONGITUDE 

«#G15.6/ 

A 

lox** longitude-azimuth 

«#G15*6) 

180 FORMAT 

! ////10X,*TERMS OF THE INCLINATION STATISTICAL EQUATIONS*/ 

A 

lOX** !I«INCLINATI0N; L-LATITUDE! Z-AZIMUTH) */ 

A 

/ lOX,* PARTIAL DERIVATIVES 

EXPECTED VALUES ♦/ 

A 

lOX** UNITLESS 


RADIANS ♦/ 

A 

10X**PIL **E15.6,* 

ELZ 

♦#E15,6 / 

A 

10X,*PIL2 **E15.6** 

EL2Z 

♦»E15.6 / 

A 

10X**PIZ **E15.6** 

ELZ2 

♦#E15.6 / 

A 

10X,*PIZ2 **E15.6** 

EL2Z2 

♦#E15.6 / 

A 

10X,*PILZ **E15.6** 

EL3Z 

♦#E15.6 / 

A 

lOX** ** 15X** 

ELZ3 

«#E15.6 ) 

190 FORMAT 

! /lOX** MOMENTS 


VARIANCE TERMS^/ 

A 

lOX,* RADIANS 


RADIANS ♦/ 

A 

10X**M3L **E15.6** 

TERMl 

♦fE15.6 / 

A 

10X**M<iL **E15.6** 

TERM2 

♦#E15,6 / 

A 

10X,*M3Z **E15.6** 

TERM3 

♦#E15.6 / 


- B5 - 



10X#«M4Z 

♦f 

E15.6,^ 

TERM4 

♦#E15.6 

/ 

lOX,* 


15X,^ 

TERMS 

♦#E15.6 

/ 

lOX,* 


15X,* 

TERM6 

♦#E15.6 

/ 

lOX,* 

♦» 

15X,^ 

TERM? 

♦»E15.6 

/ 

lOXf* 


15X»* 

TERMS 

♦#E15.6 

/ 

lOXp* 

♦# 

15X#» 

TERM9 

♦f E15.6 

) 



SUBROUTINE C ARDS ( NS AHP» NHAX ) 

C this subroutine reads fixed field data CAROS IN lA) 

C FORMAT# WRITES ON UNIT 6 FOR PRINTOUT AND WRITES ON 

C UNIT B FOR SUBSEQUENT READING IN (F) FORMAT. 

C 

IMPLICIT INTEGER(A-Z) 

DIMENSION CARD(8) 

C 

NSAMP-0 
MAX«NMAX^l 
READ(5#50) CARD 

IF (EOFC5).NE.O) GOTO 30 
URITE(6#60) CARD 
DO 10 I«1#MAX 

READ(5#50) CARO 
WRITE(6#60) CARD 
WRITE(8#50) CARO 
IF (CARD(l) .EQ.^HEND ) GOTO 20 
10 CONTINUE 

WRITE(6#70) NMAX 
GOTO AO 
20 CONTINUE 
NSAMP-I-1 

IF (NSAHP.LE.O) WR1TE(6#80) 

IF (NSAHP.LE.O) GOTO AO 
REWIND 8 
RETURN 
30 CONTINUE 
WRITE(6#90) 

AO STOP 
C 

50 FORMAT (8A10) 

60 FORMAT (5X#8A10) 

70 FORMAT ( //15X#^SUBR0UTINE CARDS - SAMPLES EXCEEDS MAX OF ♦ lA) 
80 FORMAT ( //15X#^SUBR0UTINE CAROS - NO DATA PROVIOEOA) 

90 FORMAT ( //15X#*SUBR0UTINE CAROS - NO SAMPLES INPUTS) 

END 
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c 

c 


c 


c 


c 


c 


c 

c 

c 


SUBROUTINE INPUT ( TITL El# TI TL E2 ) 

THIS SUBROUTINE READS MOOIflED NAMELIST FORMATTED DATA* 
IT READS A CARD ON UNIT 5# WRITES THE CARO ON UNIT 6# 
WRITES THE CARO ON UNIT 8 (FIRST 72 CHARACTERS ONLY)* 
THE TITLEI and TITLE2 CAROS ARE NOT WRITTEN ON UNIT 8 
BUT THE DATA IS PLACED IN THE TITLEI AND TITLF2 ARRAYS 
FOR TRANSFER BACK TO THE CALLING PROGRAM. TITLEI* AND 
TITLE2* CARDS MUST BEGIN IN COLUMN 2 WITH NO SPACES. 

THE CALLING PROGRAM MUST BLANK THE TITLE ARRAYS# 

CALL INPUTC ARG1#ARG2# )# AND RE AO ( 8# INPUTD ) • NAMELIST 
DATA MUST BEGIN WITH SINPUTD AND END WITH SEND# BOTH 
BEGINNING IN COLUMN 2. 

IMPLICIT INTEGER(A-Z) 

DIMENSION TITLE1(8)> TITLE2(8)# CARD(8) 

DATA BLANK/lOH / 


REWIND 8 
WRITE (6#40) 

10 CONTINUE 

READ(5#30) CARO 
IF (E0F(5).NE«0) STOP 
WRITE(6#60) CARD 

BLANK COLUMNS 9 AND 10 
ENCODE (10#70#W0RD) CARD ( 1 )# BLANK 
IF (WORD.NE.IOH TITLEI* ) GOTO 20 
CARO READ IS A TITLE CARD 
ENCODE (72#80#TITLE1) CARO 
GOTO 10 

20 CONTINUE 

IF (WORD.NE.IOH TITLE2* ) GOTO 30 
CARO READ IS A TITLE CARD 
ENCODE (72#80#TITLE2) CARO 
GOTO 10 

30 CONTINUE 

CARO READ IS NOT A TITLE CARO 
BLANK COLUMNS 73-80 OF DATA CARD 
ENCODE (10#90#CARO(8)) CARO(B )# BLANK 
WRITE(8#S0) CARD 

IF (CARD(l).NE.lOH SEND ) GOTO 10 


C 

REWIND 6 
CALL DATE(OAT) 

CALL TIME(TIM) 

WRITE(6#100) OAT#TIM 
RETURN 
C 

40 FORMAT (IHl ) 

50 FORMAT (8A10) 

60 FORMAT (lOX.BAlO) 

70 FORMAT (A8#A2) 

80 FORMAT (R2#7A10) 

90 FORMAT (A2#A8) 

100 FORMAT ( /////#10X#+0ATE IS *A9/ 
A 10X#+TIME IS ^49) 


END 
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SUBROUTINE C0RSI6 ( X» Y»NS AMP , SIGMAX, SIGMAY, RHOX Y, 

XMEAN»YMEAN»IMEAN» ICOR ) 

CALCULATES ME AN» STANDARD DEVIATION AND CORRELATION COEFFICIENTS 
OF TWO VARIABLES OF NSAMP SAMPLES. 

**«*««««««*** INPUT DATA ♦♦♦♦♦♦♦♦♦*♦*♦♦* 

X ■ ARRAY OF SAMPLE VALUES OF VARIABLE X 

Y • ARRAY OF SAMPLE VALUES OF VARIABLE Y 

NSAMP ■ NUMBER OF SAMPLE VALUES 
IMEAN ■ 0 IF MEAN VALUES OF VARIABLES X AND Y 
ARE TO BE SET TO ZERO. 

- 1 IF MEAN VALUES OF VARIABLES X AND Y 
ARE TO BE CALCULATED. 

OUTPUT data ♦♦♦♦♦♦♦♦♦♦•♦♦*♦ 

SIGMAX - STANDARD DEVIATION OF VARIABLE X 
SIGMAY ■ STANDARD DEVIATION OF VARIABLE Y 
RHOXY • CORRELATION COEFFICIENT OF VARIABLES X AND Y 
XMEAN ■ MEAN VALUE OF VARIABLE X 

YMEAN • MEAN VALUE OF VARIABLE Y 

IMPLICIT REAL (A-H,0-Z) 

DIMENSION X(l)fY(l) 

C 

C 

XSUM-0. 

YSUM-0. 

XSOSUM-0. 

YSOSUM-0. 

XYSUM-0. 

RHOXY-0. 

SAMPN-NSAMP 
DO 10 I«1»NSAMP 
XSUM«XSUM*X(I) 

YSUM»YSUM+Y(I) 

XSOSUM-XSOSUM«X(I)«X(I) 

YSOSUM-YSOSUM+Y(I)*Y(I) 

XYSUM-XYSUM*X(I)*Y(I) 

10 CONTINUE 
XMEAN-0. 

YMEAN-0. 

IF (IMEAN. NE.O) XMEAN-XSUM/SAMPN 
IF (IMEAN. NE.C) YMEAN-YSUM/S AMPN 
SIGMAX-SORT(XSQSUM/ (SAMPN-1. )-XMEAN**2) 

SIGMAY-SQRT( YSQSUM/( SAMPN-1. >-YME ANA*2) 

IF (ICOR.NE.O) 

A RH0XY«(XYSUM-SAMPN*XMEAN*YMEAN) /SAMPN/SIGMAX/SIGMAY 
RETURN 
END 
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SUBROUTINE STAT( DXf OYt NSAMP, ICOR# IME AN ) 

THIS SUBROUTINE CALCULATES EXPECTED VALUES AND STATISTICAL 
MOMENTS OF TWO VARIABLES. 

^^^^^t**** INPUT DATA 

OX • ARRAY OF SAMPLE VALUES OF VARIABLE X# DEG 

DY • ARRAY OF SAMPLE VALUES OF VARIABLE Y, DEG 

NSAMP- NUMBER OF SAMPLES 

ICOP • NON-ZERO VALUE RESULTS IN INCLUSION OF CORRELATION 

TERMSf I.E.# EXPECTED VALUES BETWEEN THE TWO VARIABLES 
ARE CALCULATED. 

IMEAN* NON-ZERO VALUE RESULTS IN CALCULATION OF MEAN VALUES 
OP THE TWO VARIABLES. ZERO VALUE ASSUMES ZERO MEANS. 

♦ OUTPUT DATA 

EXY - EXPECTED VALUES OF VARIABLES X AND Y 
M3X - THIRD moment OF VARIABLE X 

M3Y - THIRD MOMENT OF VARIABLE Y 

MAX - FOURTH MOMENT OF VARIABLE X 
MAY « FOURTH MOMENT OF VARIABLE Y 

IMPLICIT REAL(A-H#M#0-Z) 

COMMON /STATS/ EX Y # £X2Y#EXY2# EX3Y# EX Y3f EX2Y2 # M3 X# M3Y# MAX, MA Y 
DIMENSION DX<1),DY(1) 

DATA OEG/57. 29578/ 

C 

C 

SAMP-NSAMP 

EXY-0. 

EX2Y-0. 

EXY2-0. 

EX3Y-0. 

EXY3-0. 

EX2Y2-0. 

M3X«0. 

M3Y-0. 

MAX-0. 

MAY-0. 

XSUH-0. 

YSUM.O. 

EX-0. 

EY-0. 

IF (IMEAN.EO.O) GOTO 20 
DO 10 I-1#NSAMP 
XSUM-XSUM+DXd) 

YSUM-YSUM + DYU) 

10 CONTINUE 

EX-XSUM/SAMP 
EY-YSUM/SAMP 
20 CONTINUE 

00 30 I-l, NSAMP 

DDX-(DX(I)-EX)/DEG 

DDY-(DY(I)-EY)/OEG 

DX2-D0X«^2 

0X3-0DX^^3 

DXA-D0X*+A 

0Y2-D0Y^^2 
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OY3«ODY**3 

0Y^«0DY**4 

H3X-M3X+DX3 

f<3Y-M3Y*0Y3 

M^X«M4X*DX4 

H4Y-M4Y+DY4 

IF (ICOR.EO.O) GOTO 30 
EXY-EXY*OOX*OOY 
EX2Y«EX2Y+0X2*0DY 
EXY2»EXY2+DDX*DY2 
EX3Y-EX3Y+0X3+D0Y 
EXY3-EXY3+00X*0Y3 
EX2Y2-EX2Y2+0X2*0Y2 
30 CONTINUE 

M3X-N3X/SANP 

M3Y-H3Y/SAHP 

HAX-MAX/SAMP 

MAY-NAY/SAHP 

EXY.EXY/SANP 

EX2Y-EX2Y/SAHP 

EXY2-EXY2/SAHP 

EX3Y-EX3Y/SAMP 

EXY3-EXY3/SAMP 

EX2Y2-EX2Y2/SAMP 

RETURN 

END 
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SCIENTIFIC DATA PROCESSING ROUTINE 
SUMMARY DOCUMENTATION 


IDENTIFICATION 

Title Moments of Inclination Error Distribution (MIED) 

Routine No. 7229 Date Filed April 70 Security Class. U 

Responsible Engineer T. R. Myler 

Date Completed April 1970 Source FORTRAN Other 

Language: IV 

Key Words Orbit inclination, statistics 

RESOURCE REQUIREMENTS 

Typical CPU 10 sec Machine (s) CPC CYBER 173 No. Source Cards 519 

Core 42k (octal) Tape none Plot no Graphics none Other 

•DESCRIPTION 


Purpose : To compute inclination error statistics based on flight experience 

using a closed form solution. 

Input : Predicted insertion conditions and observed deviations of past 

flights; nominal insertion conditions of orbit requiring inclination error 
statistics. 

Output ; Mean, standard deviation, variance, third moment, skewness, fourth 
moment and kurtosis of inclination error; means, standard deviations and 
correlations of azimuth, latitude and longitude. 

Functional Description ; A Taylor series expansion is used to generate 
statistical moments of a multivariable system, which is the relationship of 
inclination to latitude and azimuth. 


DOCUMENTATION 

Vought Report 2-53030/1R-52638, "Moments of Inclination Error Distribution 
Computer Program" dated 1 March 1981. 
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1 Report No 2 Government Accession No 

NASA CR-1 65824 

3 Recipient's Catalog No 

4 Title and Subtitle 

Moments of Inclination Error 
Distribution Computer Program 

5 Report Date 

December 1981 
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T. R. Myler 

8 Performing Organization Report No 

10 Work Unit No 

9 Performing Organization Name and Address 
Vought Corporation 
P. 0. Box 225907 
Dallas, TX 75265 

11 Contract or Grant No 
NASI -15000 

13 Type of Report and Period Covered 

Contractor Report 

12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14 Sponsoring Agency Code 

490 - 02 - 02 - 77-00 

15 Supplementary Notes 

Langley Technical Monitor: R. J. Keynton 


16 Abstract This report describes a FORTRAN coded computer program which calculates 
orbital inclination error statistics using a closed-form solution. This solution uses 
a data base of trajectory errors from actual flights to predict the orbital inclination 
error statistics. The Scout flight history data base consists of orbit insertion 
errors in the trajectory parameters - altitude, velocity, flight path angle, flight 
azimuth, latitude and longitude. 

The Scout data base includes flights with launch azimuths ranging from easterly to 
southerly to slightly west of south. Latitude errors, which produce inclination 
errors, result from crossrange errors on easterly flights and range errors on southerly 
flights. Therefore, latitude errors of the flight history data base are not a 
consistent sample set of errors for a given flight azimuth at insertion. Since it is 
necessary to have a consistent set of sample latitude errors, the data base latitude 
errors are adjusted for the flight azimuth of interest. Knowing the nominal flight 
azimuth, latitude and their errors on each flight of the data base, range and cross- 
range errors are calculated. Since range and crossrange errors are independent of the 
flight azimuth, these errors can be converted to latitude errors for the flight 
azimuth of interest. 


The methods used to generate the error statistics are of general interest since they 
have other applications. Included in this report are program theory, user instructions 
output definitions, subroutine descriptions and detailed FORTRAN coding information. 


17 Key Words (Suggested by Author(s)) 18 Distribution Statement 
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FEDD Distribution 
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